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Abstract. The spectrum of the quasar PHL 1226 is known to have a strong Mgn and sub-damped Lymano- (sub-DLA) 
absorption line system with iV(Hl) = (5 ± 2) x 10' 9 cirT 2 at z = 0.1602. Using integral field spectra from the Potsdam Multi 
Aperture Spectrophotometer (PMAS) we investigate a galaxy at an impact parameter of 6'.'4 which is most probably responsible 
for the absorption lines. A fainter galaxy at a similar redshift and a slightly larger distance from the QSO is known to exist, 
but we assume that the absorption is caused by the more nearby galaxy. From optical Balmer lines we estimate an intrinsic 
reddening consistent with 0, and a moderate star formation rate of 0.5 M Q yr _I is inferred from the Ha luminosity. Using 
nebular emission line ratios we find a solar oxygen abundance 12 + log(0/H) = 8.7 + 0.1 and a solar nitrogen to oxygen 
abundance ratio log (N/O) = -1.0 ± 0.2. This abundance is larger than those of all known sub-DLA systems derived from 
analyses of metal absorption lines in quasar spectra. On the other hand, the properties are compatible with the most metal rich 
galaxies responsible for strong Mg n absorption systems. These two categories can be reconciled if we assume an abundance 
gradient similar to local galaxies. Under that assumption we predict abundances 12 + log (O/H) = 7.1 and log (N/O) = —1.9 for 
the sub-DLA cloud, which is similar to high redshift DLA and sub-DLA systems. We find evidence for a rotational velocity of 
~200 km s over a length of ~7 kpc. From the geometry and kinematics of the galaxy we estimate that the absorbing cloud 
does not belong to a rotating disk, but could originate in a rotating halo. 

Key words, galaxies: abundances - galaxies: kinematics and dynamics - galaxies: quasars: individual: PHL 1226 



| 1 . Introduction fields around QSOs which exhibit strong Lya absorption lines 

? ^ ■ have revealed several absorbing galaxy candidates close to the 

>f Metal absorption line systems in QSO spectra are fre- ,. c ■ r „ . iiVnmk c n 

n F J line of sight (e.g. Le Brun et al. 1997). Follow-up specnoscopy 

5h ' quently used to derive abundances of intervening neutral . . c ,. , . ,. , 

d ls necessary to confirm the candidates as the corresponding ab- 
. . . i gas clouds of high column density. A classic definition di- , • 

& & J sorbing galaxies, 

vides strong QSO hydrogen absorption systems into differ- 

ent regimes: Lyman limit systems (LLS) with column den- To date Qnly few {U) i o w _ redshi f t DLA s and sub-DLA 

sities of neutral hydrogen in the range 10 17 < iV(Hl) < ga i axies have been identified JTrTrnshek et alJbOOltlRao et alJ 

2 X 10 2 ° cm- 2 , and damped Lymana systems (DLAs) having fij^ Lanzett a 2003 Lacv et al]l2003l and references 
N(H i) > 2 x 10 20 cm" 2 dWolfe et al.ll986l) . Sub-DLA systems 
with 10 19 < JV(Hi) < 2 x 10 20 cm" 2 also show broad ab- 
sorption line wings characteristic for DLAs and c ontain a large 

fraction of the neutral hydrogen at high redshifts frerouxetalj | 2 QQ4. Considering the difficulties of confirming the absorb- 

|2003J>. Strong hydrogen absorption lines in QSO spectra are ing galaxieSi alternative approaches have been carried out 

believed to arise in gas-rich (proto)-galaxies, but any relation tQ determine which types of galaxies produce DLA sys- 

between absorption line systems and the galaxies responsible tems Through measurements of relative abunda nces of var- 

for these is not well understood. Deep imaging studies of the ioufj element ^ it has been suggested that DLA galaxies are 

* Based on observations obtained at the German-Spanish likely dwarf galaxies jPettini et : al Jl999[|Prochaska . et alj2003j: 

Astronomical Center, Calar Alto, operated by the Max-Planck- Pessauges-Zavadsky et alJ |200j jNissen et alJ |2004J). Sub- 

Institut fur Astronomie Heidelberg jointly with the Spanish National DLAs show the same spread in metallicity as the DLA systems 

Commission for Astronomy indi cating a similar nature of the two samples llPeroux et alJ 

Correspondence to: lchristensen@aip.de 120031) . Although metallicities of DLA and sub-DLA systems 
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can be easily derived from the QSO spectra, the metallicity 
of the parent gala xy has only been de termined for one DLA 
galaxy previously dSchulte-Ladbeck et all2004l) . 

Studies have shown that Mgn absorption systems arise 
in halos of a wide range of galaxy types ( Bergeron & Boisse 
Il99lt ISteidel et all 1 1994 and are typically associated with 
LLS. Mgn s ystems have been suggested to be related to ro- 
tating disks JCharlton & Churchilll ll998l) . while ISteidel et alJ 
(2002) found that Mgn systems cannot be explained by sim- 
ple disk models and suggested that the absorption occurs in 
rotating halos. 

The quasar PHL 1226 (Q015 1+045 at z = 0.404) has a 
strong Mgn absorption system at z — 0.1602 dBergeron et alJ 
1988, hereafter B88). A column density of N(U i) = (5 ± 2)x 
10 19 cm -2 was measured in a UV HST spectrum classifying 
the object as sub-DLA (Rao 2004, private communication). A 
bright V = 19.2+0.3 galaxy at a projected angular separa- 
tion of 6'.'4 to the west of the QSO was identified as being re- 
sponsible for the absorption system by B88, who named this 
galaxy G4. For a fiat cosmological model with Qa = 0.7 and 
Hq = 70 km s _1 Mpc~', which we use throughout the paper, 
an impact parameter of 6'.'4 corresponds to 17.6 kpc. An ad- 
ditional fainter galaxy, G3 (V = 20.5) at a similar redshift, 
was found at an impact parameter 1079 to the north of the 
QSO. Either galaxy could be responsible for the absorption sys- 
tem considering the scaling-law between Mgn halo sizes and 
galaxy luminosities ([Be rgeron & Boissg Il99lt iLe Brun et alJ 
1 19931 iGuillemin & Berg eron||l9 97|). Optical spectroscopy of 
G3 and G4 is presented in El lison et alJ ( 120041) . 

In this paper we present a study of the galaxy G4 using 
integral field spectroscopy with the P otsdam Multi Aperture 
Spectrophotometer (PMAS Roth et alJl2000l) . We describe the 
data reduction in Sect. |2] Spectra and synthetic narrow-band 
emission line images are presented in Sect. [3]From the spectra 
we derive the intrinsic reddening, oxygen and nitrogen abun- 
dances, and star formation rate in Sect.0] We analyse the kine- 
matics of the galaxy itself in Sect. [5] A discussion is presented 
in Sect.|6] and the conclusions in Sect.Q 

2. Observations and data reduction 

PMAS is an Integral Field Instrument developed at the 
Astrophysikalisches Institut Potsdam. It is currently installed 
at the Calar Alto Observatory 3.5m Telescope in Spain and 
has been available as a common user instrument since July 
1, 2002. Observations using integral field units have the ad- 
vantage of combining imaging and spectroscopy information 
in one dataset. In 2002, we started a project to use PMAS for 
the study of sub-DLA and DLA galaxies. The great advantage 
of using integral field spectroscopy for the investigation of the 
very few currently known sub-DLA and DLA galaxies is that 
they allow us to obtain data on the galaxies' positions, veloc- 
ity fields and star-forming properties all in one data cube. Such 
data are critical to derive the nature of the galaxies causing ab- 
sorption systems in QSOs. It is our goal to use this new in- 
strument to identify galaxies responsible for low- as well as 
high-redshift QSO absorbers, and to study their properties. 



During the pilot observing run for this project, we ob- 
served two QSO fields. Observations for the DLA galaxy to- 
ward Q223+131 revealed an extended Lyman-o' emission neb- 
ula surrounding the galaxy responsible for th e DLA; the analy- 
sis was presented in Christens en et al .1 (12004 . A second object 
observed in the pilot study was PHL 1226 (Q015 1+045). The 
8"x8" PMAS field was targeted at the galaxy G4 because it is 
closer to the QSO than G3, and thus the more likely absorbing 
galaxy. An additional criterion was technical - G4 is more than 
a magnitude brighter than G3 which is farther than 8" away 
from PHL 1226; its observation would have required another 
set-up and a large amount of observing time. Here, we report 
in detail on the results which can be obtained from the PMAS 
data of G4. 

We obtained 4x1 800s exposures of G4 divided over two 
nights on Sep. 5 and 8, 2002, using a grating with 300 
lines mirT 1 which resulted in a spectral resolution of 6.6 A 
FWHM. The chosen grating angle allowed to cover the wave- 
length range 4575-7880 A. The first night was photometric 
with a seeing of 1", while the seeing was varying during the 
second night and the conditions were not photometric. A log of 
the observations is presented in Tabled 

The PMAS instrument has two cameras; one for the spec- 
trograph, and an additional camera used for acquisition and 
guiding (A&G camera) which is equipped with a lkxlk SITe 
CCD. Using data from the A&G camera, one can estimate the 
variations and the evolution of the sky quality during the night. 
Photometry of the guide star images taken during the spec- 
tral integrations, show variations with a standard deviation pre- 
sented in column 5 of Tabled 

The spectrograph is coupled by 256 fibers to a 16x16 ele- 
ment micro-lens array. During the observations each lens cov- 
ered 075x075 on the sky. The detector is a 2kx4k SITe ST002A 
CCD which was read out in a 2x2 binned mode. 

Data reduction was performed with IDL bas ed routines 
written specifically for PMAS data dBeckerl Eo02). After bias 
subtraction the spectra were extracted using information of the 
location of the 256 spectra on the CCD obtained from a cali- 
bration frame obtained with an exposure of a continuum emis- 
sion lamp immediately before or after each target exposure. 
In the extraction a Gaussian line profile was used to increase 
the signal-to-noise ratio of the extracted spectra. Wavelength 
calibration was done using calibration spectra of Hg-Ne emis- 
sion line lamps also obtained following the science exposures. 
The accuracy of the wavelength calibration was checked us- 
ing sky emission lines, showing a standard deviation of 0.3 A. 
Corrections for varying fiber-to-fiber transmissions as a func- 
tion of wavelength uses flat field spectra obtained from the twi- 
light sky. 

Cosmic rays were remove d from the sp ectra using the L.A. 
Cosmic routine within IRAF ( I van Dokkuml200ll) . Each of the 
4 exposures were corrected for an average extinction value 
appropriate for Calar Alto before combining them. The ef- 
fect of differential atmospheric extinctio n was corrected us - 
ing the theoretical approach described in Filippenko (1982). 
We checked whether the data cube was appropriately corrected 
for the differential atmospheric refraction by cross-correlating 
each monochromatic image with the broad-band image shown 
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date (UT) 


exposure time (s) 


airmass 


seeing 


0- s ,ar(%) 


S/N 


Sep. 5 2002 02:48 


1800 


1.19 


1" 


0.5 


5.9 


Sep. 5 2002 03:20 


1800 


1.19 


1" 


0.5 


5.7 


Sep. 8 2002 00:03 


1800 


1.57 


1'.'3-1'.'5 


10 


3.9 


Sep. 8 2002 00:36 


1800 


1.14 


1'.'3-1'.'5 


9 


4.2 



Table 1. Log of the observations. cr star denotes the standard deviation of flux in images of the guide star during each integration. 
Column 6 lists the signal to noise ratio of a one-dimensional spectrum within a 1" radius of the QSO. 



in Fig. [2 This showed that the relative shift with wavelength 
was smaller than (X'08, which was negligible for the further 
analysis. Subtraction of the sky background was done by cre- 
ating an average sky spectrum by selecting spaxels (spatial el- 
ements) located between the QSO and the galaxy, uncontam- 
inated by flux from any of these two. This sky spectrum was 
subtracted from all 256 spectra. 

The method for combining the individual exposures was as 
follows. Firstly, monochromatic images were made from the 
data cube at some selected wavelength. Then we found a scale 
factor between the images, by calculating the total flux in the 
images. Because the spectra taken on Sep. 8 were of poorer 
quality than the first ones, the individual data cubes were scaled 
to the ones from Sep 5., and the final combination took the 
varying seeing into account by applying a weighting scheme, 
where the weights were given from the signal to noise ratio of 
each spectrum. The signal-to-noise ratio was found from a ID 
spectrum created by co-adding spectra within 1" radial aper- 
ture centered on the galaxy. These S/N ratios are listed in col- 
umn 6 in Table [0 Finally, flux calibration was done the stan- 
dard way by comparing the spectra obtained of the spectropho- 
tometric standard, BD+28°4211 observed on Sep. 5 with table 
values. After the data reduction the spectra are contained in a 
data cube of dimensions 16 x 16 x 1008 pixels. 

Further analysis of the reduce d data cube w as done us- 
ing the Euro3D Visualization Tool JSanchezl2004 . while one- 
dimensional spectra were analysed using both IRAF and our 
own software (Sanchez et al., 2004, in prep.). Fig.[0presents a 
composed broad-band image, with dimensions 8"x8". To the 
east (left), the QSO is seen at the edge of the field. The absorb- 
ing galaxy G4, identified by B88, is the principal source in the 
field at the position (-1.5, -6.5) in agreement with the impact 
parameter of 6'.'4. 

3. Galaxy and QSO spectra 

In the dust maps of lSchleeel et alJ Jl998h a Galactic reddening 
of Eb-v = 0.051 is found in the direction toward PHL 1226. 
The reduced spectra were de-reddened using this value before 
deriving other quantities. From the reduced data cube, one can 
create spectra from chosen spaxels or monochromatic- and nar- 
row band images in any wavelength region. A spectrum of the 
QSO created by co-adding 20 spaxels is shown in Fig.|3] 

One sees that strong telluric absorption lines are present 
around 7600 A which coincidentally corresponds to Ha at the 
redshift of the absorber. When line fluxes are derived, this fea- 
ture will skew the results towards lower values if not corrected 
for. Thus, to correct for this effect, a model for the telluric ab- 
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Fig. 1. Broad band image in the wavelength range 
4500 < A < 8000 A of the QSO at the position (0,0) 
and the galaxy G4 at (-1.5,-6.5) as indicated. Each spaxel 
represents a spectrum and the field of view is 8"x8". North 
is up and east is left. Coordinates are relative to the QSO 
position. 

sorption line was created. The QSO spectrum was smoothed 
using a Gaussian function with cr — 10 A, normalised to 1 at 
7580 A and outside the region around 7600 A the value was 
set to 1 . Following, all spectra were divided by this model to 
correct for the absorption as demonstrated in in Fig.|4l which 
brought out the presence of the [N n] (6548 A) line in the spec- 
trum of G4 as explained in detail in Sect. 13.21 

Co-adding all spaxels (112 in an area of 4"x7") corre- 
sponding roughly to the size of the galaxy in this dataset re- 
sults in the spectrum presented in Fig. [3] This spectrum gives 
the total continuum flux from the galaxy with S/N = 6 and does 
not show strong emission lines, specifically not before the cor- 
rection for telluric absorption is applied. One must note that 
co-adding 112 spaxels dilutes any emission line signal coming 
from a potentially smaller region. Additionally, residuals from 
the background subtraction can artificially enhance emission 
lines. This appears to be the case for the sulfur lines in Pig. El 
because they lie close to strong sky lines. In order to increase 
the signal to noise ratio of the emission line spectrum, and be 
less affected by sky subtraction residuals, a localization of the 
emission line region needs to be done. 
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Fig. 2. Spectrum of PHL 1226 created by co-adding 20 spaxels. 
The flux levels are relative since the QSO is located at the edge 
of the field of view. Broad emission lines from Hy and H/3 at 
Z = 0.404 can be seen at 6100 A and 6820 A, respectively. 
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Fig. 3. Spectrum of the galaxy created by co-addition of 112 
spaxels. Any emission line present in the galaxy has been 
smeared out because many spaxels do not contribute to the 
line emission as shown in Sect. 13. II The seemingly strong [S n] 
lines are partly caused by errors in sky subtraction. 



To check the spectrophotometry, we convolved the galaxy 
spectrum with a transmission curve of the Bessell V band filter. 
We find V - 19.08+0.05, corresponding to an absolute magni- 
tude My = -20.3 for the galaxy in agreement with the photom- 
etry in B88. 



3.1. Emission line images 

Emission line images are created by selecting appropriate nar- 
row band filters from the data cube. These narrow band images 
are created in a wavelength range of + 10A around the emis- 
sion wavelength as shown in Fig. [5] The continuum emission is 
subtracted using narrow -band images in wavelengths adjacent 
to the emission lines. Compared to the broad band image in 
the lower right panel, one sees that most of the emission lines 
originate ~1" to the south of the G4 centre. Furthermore, the 
strongest emission is located in a region within a 2"x2" aper- 
ture. No other strong emission line regions are found in the 
field. 




7400 7500 7600 .7700 

Wavelength [A] 



7800 



Fig. 4. Zoom in of the spectrum in Fig. [6] in the region around 
[Nn], Ha, and [S n]. This plot shows the effect of correcting for 
the telluric absorption feature around 7600 A. The lower spec- 
trum (thin line) is the uncorrected spectrum, and the gray line 
is the model of the telluric feature created from the QSO spec- 
trum. The upper spectrum (thick line) has been corrected for the 
absorption feature and offset by +3xl0~ 17 ergcrrT 2 s A -1 
for clarity. 

3.2. High surface brightness region 

To create a spectrum with a higher signal to noise ratio relevant 
for deriving properties from the emission lines, we selected 
spaxels from Fig. [5] where the surface brightness of the emis- 
sion lines is strongest, i.e. within a 1" radial aperture (corre- 
sponding to ~15 spaxels) centered on (-2.5,-6.5). Several lines 
can be identified in the spectrum shown in Fig. [6] Observed 
line fluxes and equivalent widths (EW) are listed in Table [2] 
All emission line fluxes refer to the spectrum which has been 
corrected for the telluric absorption feature. 

The median redshift of the emission lines is z e m = 0.1595 + 
0.0006 in agreement with z = 0.1592 + 0.0020 found by B88. 
The observed FWHM of emission lines are 7.2 + 0.5A and 
6.1 + 1.1 A for Ha and H/3, respectively, which implies that 
the lines are consistent with being unresolved given the spec- 
tral resolution. Line fluxes of the faint lines have been derived 
assuming a similar width as for the brighter lines. The velocity 
difference between the absorption system and the median emis- 
sion redshift is 180 km s _1 . In Sect. |5] we will return to further 
analysis of the kinematics involved in this system along with 
the implications. 

To check if the correction for telluric absorption is appro- 
priate we measure the line r atio [Niil/t 6583/[Nnl/t6548 which 
should theoretically be 3.0 JWiese et al.ll 19661) . while the de- 
rived ratio is 3.58 + 0.76. Thus, we conclude that the correction 
does not introduce large errors of the derived line fluxes. Before 
application of the correction the line ratio is ~2, however, it is 
uncertain as the continuum level is wrongly placed because of 
the telluric absorption. 

Although unaffected by the correction for telluric absorp- 
tion, the [Om]/*5007/[Om]/14959 line ratio should be 2.87, 
while we measure the ratio 2.31+0.85. Line ratios of the emis- 
sion lines [Om]/I5007 vs. U/3 and [Nn]A6583 vs. Her clas- 
sify this o bject as an H n galaxy accord ing to the classification 
scheme of Veilleux & Osterbrock ( 1987). 
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Fig. 5. Emission line images of various emission lines where the continuum has been subtracted using images using adjacent 
wavelengths. For comparison a broad band image is included in the lower right hand panel. Smoothed contours of 2, 4, 6,...cr 
levels are overlayed. The solid black line represent the 4<x level. 



Using long slit spectroscopy B88 find a Hfi line flux of 
2.2 x 1CT 16 erg cmT 2 s _1 , which is slightly below our measured 
line flux given in Table|2] Their long slit spectrum was obtained 
with a l'/5 slit, i.e. smaller than our 1" radial aperture for cre- 
ating the spectrum. Therefore, our finding of a larger line flux 
can be explained by small slit-losses in their spectrum. 

4. Derived properties 

Using the strong emission lines observed in the galaxy spec- 
trum we derive the dust reddening, abundances from the nebu- 
lar emission lines in Table|2 and its star formation rate. 

4.1. Dust reddening 

From the observed Balmer decrement the dust content can be 
estimated. In the Case B Balmer recombination scenario the 
ratio of the emission lines is I(Ha)/I(H/3) = 2.85, d epending 
slightly on the temperature and density dBrocklehurstlll97ll) . 
while the observed flux ratio is F(Ha)/F(Hfl) = 3.02+0.39 
which is consistent with no reddening. 

Us ing the Milky Way extinction curve in iFitzpatrickl 
Jl999l) . we find the reddening E B - V = 0.05+0.09 or 
equivalently Ay = 0.16+0.28 for a Galactic value of 
R v = 3.1. Similarly, the theoretical ratio of I(Uy)/I(Ufi) 
is 0.46 in the case B scenario, while we measure 
F(Hy)/F(Hf3) = 0.52 ± 0.17. This corresponds to 



Eb-v-0- 25+0.88 which is in agreement with the reddening de- 
rived from the HafH/3 line ratio. 

Because the internal extinction in this system is consistent 
with 0, we do not correct the derived emission line fluxes for 
this effect. With integral field data one can in principle create 
dust maps by dividing the Ha-image with the Fl^S-image, but in 
this case, the signal to noise ratio is not high enough to derive 
a reliable extinction map. 

4.2. Chemical abundances 

Using the line fluxes of strong emission lines we derive abun- 
dances of oxygen and nitrogen using various calibrations and 
diagnostics from the literature. 

4.2.1 . Oxygen 

To determine the oxygen abundance of the 
galaxy we make use of the relation 03N2 = 
log qO lliU5007/Hg) /(fN IIU6583 /Hor) recently calibrated 
in IPettini & Pagell J2004I) . Our data cover all the strong 
lines involved in the 03N2 calibration which for G4 yields 
12 + log(0/H) = 8.66 + 0.10, while an additional error of 0.14 
dex is due to the scatter in the calibration itself. The abundance 
corresponds to the solar oxygen abundance 8.66+0.05 in 
Aspl und et alJ (120041) . This solar oxygen abundance is lower 
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Fig. 6. Spectrum of the galaxy created by co-adding 15 spaxels that appeared to be associated with the bright line emission region. 
Emission lines listed in Tableware indicated. 



line 


A" 


redshift 


line flux 


EW 




(A) 




xl0~ 16 erg cirr 2 s _I 


(A) 


Hy (4340.49) 


5037.29 


0.1602 


1.53+0.48 


3.8 


BJ3 (4861.36) 


5639.60 


0.1598 


2.97+0.37 


8.9 


[Om] (4958.91) 


5754.00 


0.1600 


0.70±0.23 


2.8 


[Om] (5006.84) 


5808.78 


0.1598 


1.62±0.28 


6.6 


[Nn] (6548.05) 


7594.64 


0.1595 


0.83±0.17 


7.6 


Ha (6562.85) 


7611.88 


0.1595 


8.98±0.28 


43.8 


[Nn] (6583.45) 


7635.31 


0.1595 


2.97±0.17 


16.2 


[Sn] (6716.44) 


7788.96 


0.1594 


3.04±0.60 


17.1 


[S n] (6730.82) 


7804.33 


0.1592 


1.95+0.50 


15.5 



Table 2. List of emission lines in the galaxy spectrum and the corresponding redshifts. Line fluxes are derived after a correction 
for Galactic extinction and telluric absorption was applied. Fluxes and EWs are the observed ones. " Heliocentric vacuum values. 



than found in models previously , driving the derived oxygen 
abundance for G4 relative to solar towards higher values 
than reported in the literature for other Mg n galaxies. Taking 
inste ad th e solar oxygen abundance 8.74+0.08 obtained by 
iHolwegeil 1200 lb . G4 has sub-solar metallicity (O.8Z ), but 
still consistent with solar within lcr errors. 

The 03N2 ratio benefits from the fact that the involved 
lines are not separated by long wavelength ranges, and thus 
the quantity is largely unaffected by dust obscuration. At any 
rate, the small intrinsic reddening inferred for G4 implies that 
the 03N2 ratio is affected little. 



1 The large r solar O abundance in lHolwegeil 1200 lh compared to 
that derived in Asclund et al. 1 2004) is due to an ignored contribution 
from a Ni blend, but also a difference between their adopted ID and 
3D models. 



For comparison we also calculate oxygen abundances 
using line diagnostics calibrated by other authors. For ex- 
ample, the oxygen abundance can be estimated from the 
ratio N2 = log( r N III /Hot) using the calibration in 
iDenicolo et al.l d2002l) . These two strong emission lines 
are present in the G4 spectrum, and the calibration gives 
12 + log(0/H) = 8.77 ± 0.06. Appl ying instead the slightly 
different calibration of the N2 ratio in IPettini & Pagell (|2004), 
we find 12 + log(0/H) = 8.65 + 0.03 in agreement with 
the 03N2 diagnostic. The intrinsic lcr scatter of this latter N2 
calibration causes an additional uncertainty of 0.18 dex. 

As another line of enquiry, we combine our data of 
G4 with those available in the literature. Oxygen abun- 
dances have frequently been determined using the R23 = 
(TO Hi/13727 + JO IUU 4959 + [Om]/l5007)/H/3 intensity ra- 
tio IPa gel et alJ Il979l) . Although the spectrum of G4 does 
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Diagnostic 12 + log(0/H) 

03N2 (1) 8.66 ±0.10 

N2 (1) 8.65 ± 0.03 

N2 (2) 8.77 ± 0.06 

R 23 (3) 9.02 + 0.13 

Table 3. Oxygen abundance determin ations. Diagnostics hav e 
been taken from followin g p apers: (DlPettini &Pagell (|2004), 
(2) lDenicol6 et alJd2002l) . (3) lKobulnickv et al.N l999) 



not cover the wavelength of [On], we can estimate roughly 
the value if it had been covered. If the flux reported in B88 
([On] = 4.1 x 10~ 16 erg ctrT 2 s~') is not corrected for 
Galactic extinction, the un-absorbed value will be / D b s = 4. 1 X 
10 o.4£ B _„«(5963A)^ where Eb v is the Galactic extinction, and 
7?(5945A) = 4.125 is the value of the extinction curve at the 
wavelength of [On] at z = 0.1595. Furthermore, we correct 
the slit-loss present in their spectra by applying a scale fac- 
tor between their H/3 line flux and ours. We estimate the [On] 
line flux F(\Oii\) = (6.7 + 0.4) x 10~ 16 erg cm -2 s -1 . 
Ellis on et all 1*2004 ) find an [O n] flux twice this value, but also 
their H/3 and [Oin] fluxes are twice the ones we measure. 

Using the scaled B88 [On] line flux we find 
log Rn _ 0.48+0.18, which according to the calibrations in 
iKobulnickv etaD J 1999b gives 12 + log(0/H) = 9.02 + 0.13. 
This value agrees within the errors with t he 7 ?7vbased 
12 + log(0/H) = 8.9 + 0.2 in lEllison et alJ d2004l) which 
justifies our assumption of applying a scaling factor of the 
[On] flux. 

All the derived abund ances using the various diagnostics 
are summarised in Table|3] Kobulnickv et al. ( 1999) report that 
the calibration is very uncertain for line fluxes with lower sig- 
nal to noise ratio than 8:1. Given such uncertainties of the lat- 
ter abundance determination, combined with the fact that the 
03N2 calibration involves emission lines detected in the data, 
and not a scaled [On] flux derived from B88, we rely on the 
oxygen abundance derived from the 03N2 ratio. Thus, we find 
that the galaxy has a solar oxygen abundance. 

4.2.2. Nitrogen 

Having derived a solar oxygen abundance for G4, this in- 
formation can be used to estimate the electron temperature 
which in turn is used to determine the N abundance. Because 
oxygen is the main coolant of a gaseous nebula, one ex- 
pects to see a correlation between the oxygen abundance and 
the electron temperature of an Hn region. From a sample 
of ext ragalactic Hn regions with a large range of O abun- 
dances fvan Zee et all (119 98b) derive T e . In their 39 regions with 
8.56<12+log(0/H)<8.76, similar to the value of G4 within er- 
rors, we find an average T e = 7000 + 900 K. This temperature 
estimate is in agreement within lcr err ors with the rela t ion be - 
tween R23 and the [Nn] temperature in lThurston et alJ (TT996). 
which yields f[Nn] = 7600 + 1000 K. The uncertainty both re- 
flects the calibration scatter (500 K) and the uncertainties of 
line fluxes (900 K). 



Diagnostic log(N/0) 

"(1) -0.98 ± 0.23 

(2) -0.90 ± 0.22 

(3) -0.98 ±0.13 

(4) -0.90 ± 0.08 

Table 4. Nitrogen abundance determination s using different di - 
agnostics f r om fol lowin g papers: (1) Ivan Zee et alJ dl99 8b): 
IPagel et alJ dl992l). (2) iThurston et alJ dl996b: IPageletal 
|l992h. (3) IShaw & Dufourl ill 9951) . (4) iKewlev & Donite 



Abundances of ionized nitrogen can then be derived using 
the relation betwe en temperature and nitrogen to oxygen ratio 
(Pagel et al. Il992l) . Disregarding ionization corrections, the ni- 
trogen abundance can be derived assuming (N/O) = (N + /0 + ), 
which is a valid approximation since the ionization potentials 
for Q + and N 4 " are similar. Fu rthermore, as shown by the mod- 
els in iThurston et alJ Jl996l) . this approximation is accurate 
within 5%. For the emission lines from G4 we find log (N/O) = 
-0.98+0.23 using the average T e , i.e. sub-solar, but consistent 
within lcr with the solar value log (N/O) = -0.81 dHolwegerl 
200 if). A similar result is obtained from the [Nn] tempera- 
ture which gives log (N/O) = -0.90+0.22. A linear relation 
between the O and N abundance at hig h oxygen abundance s 
is observed in extragalactic Hn regions dvan Zee et alJ[l9 98a). 
thus a solar N/O abundance is expected for G4. 

Ionic abundances can also be derived using software for 
analysis of e mission line nebulae w ithin the IRAF/STSDAS 
environment dShaw&DufouJll995l) . As inputs we use the 
[N n] temperature estimated above and a low density environ- 
ment (n = 10 cm 4 ) which is preferred from the observed 
[Sn] /1/16717,6731 line ratio. These values give log (N/O) = 
-0.98+0.13 where the error mostly depends on the uncertainty 
of the temperature. Choosing instead a density of 100 cm -3 
only increases log (N/O) by 0.01 dex. 

As a consistency check, we investigated whether the de- 
rived oxygen abundance is in agreement with calibrations to 
derive the (N/O) abundance ratio. Using the polynomial re- 
latio ns between the ox ygen abundance and log(N + /0 + ) in 
iKewlev & Dopital d2002l) yields log (N/O) = -0.90 + 0.08 in- 
ferred from the oxygen abundance derived from the 03N2 
diagnostics. The error includes the uncertai nty of the chosen 
ionization parameter. However, as shown in Kewlev & Dopita 
(2002) this diagnostic is relatively independent of the ioniza- 
tion parameter for metallicities larger than half solar, which is 
the case here. The derived abundance ratios using the various 
estimators are summarised in Table |4] From these values we 
conclude that the galaxy has log (N/O) — -1.0+0.2. 

4.3. Star formation rate 

The derived line fluxes of emission lines can be used to de- 
rive an overall star formation rate (SFR) of the galaxy. A 
redshift of z = 0.1595 corresponds to a luminosity dis- 
tance of 2.35xl0 27 cm (760 Mpc) for the given cosmologi- 
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cal model. Thus, the Ho- line flux corresponds to a luminosity 
L(Ha) = (6.25 + 0.20) x 10 40 erg s -1 . Using the relation 

SFR = 7.9x 10~ 42 L(Ha) (ergs" 1 ) (1) 

in iKennicuttl Jl998l) to convert the flux to a SFR we find 
SFR = 0.49+0.15 M yr _1 . The uncertainty includes the one 
from the line flux and a larger additional uncertainty from the 
intrinsic scatter of the calibration of the conversion factor of 
-30%. 

The Ha based SFR relies only on the data set of G4 pre- 
sented here, while an alternative measure of the SFR can be 
estimated from the scaled [On] line flux. At z — 0.1595 the lu- 
minosity is L([0 III) = (4.67 +0.28 ) x 10 40 erg s _I , which, using 
the conversion in Kennicutt ( 1998), yields SFR = 0.43 + 0.12 
M Q yr -1 . This value is f urthermore in agreement with the cali- 
bration in lKewlev et ai1ll2004l) . which yields SFR = 0.45+0.17 
M yr~ 1 . The agreement between the SFR of the galaxy derived 
from different calibrations gives credibility to the inferred small 
internal reddening in the galaxy. 

5. Kinematics 

From the redshift difference between the sub-DLA cloud and 
the galaxy found from the one-dimensional spectrum of G4, 
we find that the relative velocity difference is 180 km s . 
Considering that the galaxy seems to be a fairly luminous 
galaxy, such a velocity difference over a distance of 17.6 kpc 
could be consistent with rotation of a massive disk, but the ori- 
entation of the galaxy with the major axis oriented roughly 
north-south, while the sub-DLA cloud lies to the north-east, 
does not support this hypothesis. We analyse here only radial 
velocities but note that the proper motion of the cloud could be 
significant. 

By fitting the strongest Ha emission lines with Gaussians 
we examine the dynamics of the galaxy G4. Because each in- 
dividual spectrum is rather noisy, the fitting could only be done 
satisfactorily in the region around the strongest emission. In 
Fig-El velocity offsets relative to the systemic z - 0.1595 are 
shown as contours overlayed an image of the Ha emission in- 
tensity. Uncertainties of the centroids of the Gaussian fits to the 
strongest lines are of the order of one tenth of the spectral reso- 
lution, i.e. 0.6 A corresponding to errors of 25 km s _1 . We find 
evidence for a systematic velocity of ~ -120 to +160 km s _1 
with the rotation axis oriented roughly east-west. The point of 
zero-level velocity is located roughly 01' 5 to the south of the 
centre of the G4 continuum emission. 

In this representation the absorbing cloud towards the QSO 
(at coordinates 0,0) has a velocity offset of +180 km s with 
respect to zero velocity of the Ha emission. We conclude that 
it is unlikely that the absorbing cloud is participating in a disk 
rotation. 

In the integral field data we measure a velocity of 
160 km s _1 relative to the zero velocity contour of Ha at a dis- 
tance of ~2" corresponding to 5.5 kpc. However, if one takes 
the centre of the galaxy as the true reference point, the velocity 
difference is ~200 km s _1 over a distance of 2'.'5 (6.9 kpc) as in- 
dicated in Fig. [7] Assuming that the galaxy is viewed edge-on, 
the dynamical mass is Md yn = v c r/G = (6±2)x 10 10 M G , where 



o 




ARA (arcsec) 



Fig. 7. The gray scale image shows the intensity of the Ha 
emission line. Relative velocities in km s _1 inferred from Ha 
emission lines are shown as contours. Grey contours (orange in 
the electronic edition) are positive velocities, and black (blue) 
contours are negative. Only regions showing Ha emission have 
been included for clarity. In comparison, the relative velocity 
difference of absorption in the sub-DLA cloud is +180 km s _1 
towards the QSO sight-line. 

the error is due to the uncertainty of the velocity. Because the 
inclination axis is not known, this value is a lower limit, while 
another uncertainty is due to the fact that we probably do not 
reach asymptotic values for the rotation curve. In any case, the 
estimate shows that the galaxy is a fairly massive one. 

6. Discussion 

Having derived various properties of the sub-DLA galaxy, we 
proceed by comparing them to those of other sub-DLA and 
DLA galaxies in the literature. We emphasize once again that 
there is another galaxy, G3, at the same redshift as G4 but hav- 
ing a larger impact parameter. G3 might contribute to the ab- 
sorption seen in PHL 1226. No galaxy in Fig. 3 of B88, how- 
ever, is closer to the QSO than G4, and we discuss in section 
6.3 that our data allow us to exclude an emission-line galaxy 
with a smaller impact parameter than G4 between G4 and the 
QSO. As is common practice in the literature on DLA galax- 
ies, we proceed with the assumption that the absorbing cloud 
belongs to the galaxy with the smallest impact parameter. 

6.1. Relationship with DLAs, sub-DLAs and Mg n 
systems 

The derived oxygen abundance of the galaxy G4 is larger 
than abundances derived for h ig h redshift s ub -DLA s 
JPessauges-Zavadskv et"ai1 120031 IPeroux et ail 120031) . 
Furthermore, compared to the nearby DLA absorber 



L. Christensen et al.: The candidate sub-DLA galaxy toward PHL 1226 



9 



SBS 1543+594 dSchulte-Ladbeck et alJ Eool . this galaxy 
has a higher abundance, but that goes in hand with its higher 
luminosity. In fact we calculate M B = -20.0 using the 
spectrum in Fig. [3] shifted to the rest frame and convolved 
with the Bessell B band transmission function. G4 exhibits 
values which are ent irely consistent with the local luminosity- 
metal licity relation ( Kobul nickv & Kool l2000t IPilvuein et alJ 
12004 . 

On the other hand, G4 has an oxygen abundance com- 
patible with the upper end of abundances derived for the 
Mg n selected galaxies (-0 .6 < [O/H] < -0.1 derived in 
iGuillemin & Bergeron! 19971) . For comparison, high resolution 
spectra of strong Mg n absorption systems towards a couple of 
QSOs at z ~ 1 have indicated metallicities ranging from 0. 1 to 
1 times solar (Din g et al.ll2003albr) . This could indicate a gen- 
eral agreement between metallicities derived using absorption 
lines and emission line diagnostics. 

The kinematic investigation showed that the absorber is 
clearly not associated with a disk, but it could be part of a ro- 
tating halo as suggested for other Mg n absorbers ( S teidel et alJ 
120021) . The galaxy G4 and its sub-DLA system could thus be 
similar to the z = .16377 sub-DLA towards Q0850+4400 
jLanzetta et aljfl997i) . They showed that the absorption arises 
at large galactocentric distance (17 kpc) and does not partici- 
pate in a general disk rotation. 

Alternatively. iBond et alJ J200ll) suggested that a substan- 
tial fraction of very strong Mg 11 systems arise in superwinds 
from galaxies. Along similar lines, expanding superbubbles 
were suggested to be the explanation for metal absorption line 
mirror sym metries seen in the strong Mgn absorber towards 
Q1331 + 17 fellison et alJl2003l) . A superwind from the galaxy 
G4 is unlikely at present time given the relatively small SFR. 
However, as it is comparable to an L* galaxy, a previous star- 
burst could have expelled neutral gas clouds enriched by met- 
als. 

6.2. Abundance gradient effect 

Metallicity studies of DLAs and sub-DLAs are trying to ex- 
plain the question of what is the nature of the parent galaxy. 
Only in very few currently known cases, where the parent 
galaxy has been clearly identified, can such an investigation 
be carried out. Abundances derived from metal absorption 
lines associated with the DLAs can be compared to abun- 
dances based on emission line diagnostics for Hn regions in 
the galaxy. The DLA galaxy SBS 1543+594 is one example, 
where the impact parameter is small - in fact, the sight line to 
the QSO goes through the disk of the dwarf galaxy. In such a 
case, the inferred metallicity of the DLA cloud and the galaxy 
is expected to be similar. However, in the case of the G4/sub- 
DLA system towards PHL 1226 where the impact parameter 
is larger and the sub-DLA cloud possibly does not belong to 
the disk of the galaxy, one would expect to find lower (O/H) 
and (N/O) abundances of the sub-DLA cloud due to abundance 
gradients. 

Analysing abundan ces in extragalacti c Hn regions 

IVila-Costas & Edmunds! i 19921) and |Zaritsky e t aD 1 19941) 



found a large scatter in the abundance gradients for individual 
galaxies of a given Hubble type. They also found a tendency 
for more shallow gradients, expressed in dex/kpc, for early 
and late type spirals compared to intermediate type spirals. 
Locally, three face-on disk galaxies were found to have strong 
abundance grad ients along the major axes in the outer regions 
of their disks ( Fergus on et alJ fl998. hereafter F98). Fewer 
studies have analysed the metallicity gradient along the minor 
axes of galaxies, but smaller metallicities have been found for 
extrapla nar Hn regions comp ared to the core of an edge-on 
galaxy ( Tullma nn et alJ 12003). The study of Hn regions in 
face-on galaxies can therefore also be biased by projection 
effects. 

The data presented here do not allow us to estimate the 
Hubble type of G4, making it difficult to predict an appropriate 
abundance gradient. If one takes the measured metallicity gra- 
dient from F98 with an average in log(0/H) of -0.09 dex kpc 1 , 
the abundance of the sub-DLA cloud at 17.6 kpc is expected 
to be ~1.6 dex lower than what we find for the G4 disk. 
Similarly, the average gradient in log(N/0) is -0.05 dex kpc" 4 
which implies an abundance ratio smaller by 0.9 dex. Thus 
one could expect 12+log(0/H) = 7.1 and log(N/0) = -1.9 
for the sub-DLA abundances. We assume that the gradients 
are s traight lines which may not be the case JZaritskv et alJ 
Il994l) . If a low abundance of the PHL 1226 sub-DLA cloud 
should be confirmed by future space based spectroscopy, these 
values would place the sub-DLA system toward PHL 1226 
among the metallicities for the currently measured high red- 
shift DL A and sub-DLA systems measured by several au thors 
JPettini et all2Q02bULot)ez et alJ20oilLopez & EllisorJ2003l: 
Centurio n et all2003l) 

These approximate cloud metallicities are crude estimates 
since the impact parameter is a lower limit due to the unknown 
inclination angle, and the individual abundance gradients in 
F98 vary within a factor of 2. Furthermore, the lowest metallic- 
ities observed by F98 reach 12+log(0/H) = 7.95, thus we are 
extrapolating their metallicity gradient. 

6.3. Sub-DLA cloud properties 

Other properties of G4 such as HB luminosity, EW of HB , ab- 
solute magnitude, and oxygen abundance are compatible with 
those of emission line field g alaxies at redshifts 0.26 < z < 
0.82 (Kobulni ckv et alj|2003l) . This apparently seems to be in 
contradiction with the spectroscopic analyses suggesting that 
sub-DLAs/DLAs are chemically less evolved than star form- 
ing galaxies at similar redshifts. Specifically, sub-DLAs/DLAs 
have sub-solar metalli city and element abundances suggest- 
ing low SFRs JPettini et alJ l2002at IProchaska et alJ l2003t 
Pero ux et al J 12003: Dessauges-Zavad skv et alJ 12003). On the 
other hand we can not, with the currently available data set, ex- 
clude the possibility that a galaxy less luminous than G4 closer 
to the line of sight towards PHL 1226 is responsible for the sub- 
DLA. In the integral field data presented here we find no evi- 
dence for line emission closer to the line of sight than 5" to the 
west of the QSO. With the presented single pointing towards 
the QSO we can not say anything about the other directions. 
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As argued, either galaxy G4 or G3 could be responsible 
for the sub-DLA cloud. If the cloud were associated with G3, 
it too has to be associated with a halo because of the orien- 
tation of the galaxy which suggests an elongation in the east- 
west direction (i.e. perpendicular to the direction toward PHL 
1226). Assuming that the sub-DLA cloud indeed belongs to the 
galaxy G4, we find a velocity difference of 180 km s _1 from 
the sub-DLA cloud redshift to the Ha velocity at the centre of 
G4, for which we estimate the galaxy mass (6 + 2) x 10 10 M G . 
In this case, the escape velocity at a distance of 17.6 kpc is 
200+30 km s implying that the cloud could be gravitation- 
ally bound. 

Yet another possibility for the location of the sub-DLA 
cloud is gravitational interaction between systems. Indeed the 
galaxies G4 and G3 have a distance of 12'.'6 from each other 
corresponding to 35 kpc, so it could be an interacting system, 
but with the currently available observations we can not test 
this scenario. 

7. Conclusions 

Using integral field spectroscopy with PMAS we have ob- 
serve d the absorber G4, previously identified by Berg eron et alJ 
dl988l) . toward the QSO PHL 1226 at an impact parameter of 
6'.'4. This galaxy is most probably responsible for a strong Mg n 
and sub-DLA absorption system at z — 0.1602 in the QSO 
spectrum. At the same redshift another galaxy, G3, is present. 
We cannot determine whether the PHL 1226 absorber belongs 
to galaxy G4 or G3, but concentrate on G4. Their impact pa- 
rameters of 6'.'4 and 10 / /9 correspond to 18 and 29 kpc, respec- 
tively, implying that either could be responsible for the absorp- 
tion given the scaling-law of Mg n absorbers. 

A strong emission line region is shown to be limited to an 
area of approximately 1" in radius within the galaxy. In the 
spectra we find emission lines from [Om], [Nn], [Sn] as well 
as Balmer lines Hy, U/3 and Ha at the redshift 0.1595+0.0006. 
We do not find regions of emission at the same redshift closer 
to the QSO line of sight. 

From the Balmer line ratios we find evidence of an intrinsic 
reddening of Eb-v = 0.05 + 0.16, i.e. consistent with 0. From 
the measured Ha line flux we derive a SFR=0.5 M R yr' 1 . 

Us ing the 03N2 line ratio diagnostics from lPettini & Pagel 
(2004J) we derive a solar oxygen abundance 12 + log(0/H) = 
8.7 + 0.1. Also a solar value of the abundance ratio log(N/0) 
= -1.0+0.2 is found implying a metallicity comparable to the 
upper-end of metallicities for the currently known sample of 
Mgn galaxies. 

A kinematic analysis of the Ha emission line showed that 
the galaxy has rotational velocities of -120 to +160 km s _1 
relative to the systemic redshift with the rotational axis oriented 
roughly east-west. The sub-DLA cloud, on the other hand, has 
a velocity difference of 180 km s _1 relative to the galaxy and 
an impact parameter of 17.6 kpc above the disk assuming that 
the disk is seen edge on. With such geometry and kinematics, 
the sub-DLA cloud is likely part of a rotating halo and possibly 
gravitationally bound. From relative velocity measurements we 
derive a kinematic mass of 6xl0 10 M , which corresponds to a 
fairly massive galaxy. The absolute magnitude is M v = -20.3 



and Mb = -20.0 which is consistent with the mass-luminosity 
relation for spiral galaxies (Forbes 1992). 

With future UV space-based spectroscopy it will be possi- 
ble to compare the metallicity of the sub-DLA cloud towards 
PHL 1226 with abundances derived for the galaxy, which is 
necessary in order to understand the relation between the ab- 
sorption lines in QSO spectra and the galaxies responsible for 
them. If there is a difference in metallicity in line with the 
metallicity gradients observed in local disk galaxies, we ex- 
pect that the properties of the sub-DLA cloud will be similar to 
those of high redshift DLA and sub-DLA systems. This is an 
intriguing prospect which could suggest that the specific sight 
line through the galaxy responsible for the DLA or sub-DLA 
cloud has important consequences on the derived properties of 
the cloud. 
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